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INTRODUCTION 


This Physics Lab Book was developed to provide you the opportunity to use a 
computer as a problem solving tool. You will write computer programs which will 
enable you to investigate the physics field of geometric optics. Using just one pro- 
gram, you will be able to perform many different experiments and make your 
own generalizations. 





To use Geometrical Optics, you will need the following. First, you should have a 
background in algebra and some trigonometry. Ordinarily the material covered in 
Geometrical Optics requires differential calculus, but since you will have a com- 
puter to do your calculations, you will not need this. Secondly, the Lab Book 
assumes that you already know how to write a simple computer program using 
the BASIC language. If you do not already have some experience in programming 
using BASIC, you will want to study this before you attempt the material. 
Consult the BASIC Manual for the computer you are using. Last, use of this Lab 
Book requires that you have access to a computer for at least two hours per week. 
If more time is available, you may be able to experiment further on your own, 
either to improve your program or to investigate other aspects of physics that 
interest you. 


As you will discover, there is no one ‘“‘right”’ way to use a computer as a problem 
solving tool. There are many different ways to solve one problem by program- 
ming. Experiment and learn as you go. You’ll find you are learning something new 
each time, both about your subject matter and about using the computer to solve 
problems. 


This book was designed to help you by providing several different kinds of 
material. First there are the exercises with the preparatory explanatory material 
and the drawings. These exercises are sequenced so that you can apply what you 
have learned in the previous problem in solving the next one. Often you can take 
your preceding program and simply add to it to create a program that will pro- 
vide answers to the more general or more advanced problem. 


An example program and flow chart follow the regular exercises. You may wish 
to review this flow chart and program before you begin using the Lab Book. The 
example demonstrates planning a solution (the flow chart) and the coding of the 
solution (the program). When you do begin using the Lab Book, you may choose 
to flow chart your solutions first. This is good programming practice. Drawing the 
flow chart provides a check of your logic, and the finished flow chart can then be 
used as a guide in each step of writing your program. 


The section of advanced problems is provided for any students interested in 
further work in this area. You may wish to tackle this after you have completed 
Exercise 20. 
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GEOMETRICAL OPTICS 


Why can you see around a corner using a mirror? Why does a pencil standing in a 
glass of water appear to “‘bend”’ at the surface of the water? Why does a magni- 
fying glass concentrate the light from the sun into an intense spot that is hot 
enough to ignite paper? We have all seen these things at one time or another and 
probably have wondered how or why they happen. 





In general, physics is concerned with the explanation of occurrences such as these 
in terms of general laws or ideas. The better the explanation, the more things that 
can be explained. The examples above all fall within the part of physics known as 
geometrical optics. Geometrical optics is founded upon the twin ideas of reflection 
and refraction. If we can find a satisfactory explanation of reflection and refrac- 
tion we will have a theoretical basis to understand most of geometrical optics. 


There is an idea that can be used to discover the laws of reflection and refraction. 
It turns out that variations on the same idea can be used to explain many topics 
in mechanics, another branch of physics. This idea is generally not discussed in 
introductory physics texts. The reason is because the usual mathematical treat- 
ment of the idea requires advanced techniques that are past the capabilities of the 
beginner. However, with the computer we can perform an end run around the 
mathematical difficulties and examine the idea and its consequences. 


Before proceeding, we should review some concepts about coordinate geometry. 
If you do not need this review, go on to the next section. 


GEOMETRICAL OPTICS 
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We can locate a point on a coordinate system in terms of an ordered pair of 
numbers (x, y). In the diagram in Figure 1, the point A is located by x = 2 and 
y = 1. Likewise, B is located by x = 5 and y = 5. On such a coordinate system 
the values of x and y may be either positive or negative. 





Figure | 


GEOMETRICAL OPTICS 
Coordinate Geometry 








PHYSICS 
Hewlett-Packard Computer Curriculum 


Exercise 1 — Coordinates 


In the spaces provided in the table, record the coordinates of the points in the 
diagram in Figure 2. 





Figure 2 


10 


Suppose we wanted to find the distance between two points. In Figure 3 the 
point A has coordinates (3, 2) while B is at (6, 6). We can construct a triangle as 
shown to find the distance between A and B. We use the familiar Pythagorean 
Theorem which states that the square of the hypotenuse of a right triangle is 
equal to the sum of the squares of the other two sides of the triangle, or 


(AB)* = (AC)? + (CBY? 





Figure 3 


We can find the desired length AB by taking the square root of both sides to get 
AB = V (AC)? + (CB)? 


In this example AC = (6 — 3) = 3, CB = (6 — 2) = 4. If these values are substituted 
into the equation for AB we find 


AB=V (3)? + (4% =V9 + (62x 1525 


GEOMETRICAL OPTICS 
Coordinate Geometry 
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But AC is the difference in x coordinates of the two points, while CB is the differ- 
ence between y coordinates. If we locate one point with the coordinates (x,, y1), 
and a second point at (x2, y2 ), you can easily see that the distance d between the 
two points is 


d=VJ(x2 - x1)? + (y2 -y1) 


Exercise 2 — Distance Between Points 


Write a program which will compute the distance between two points. The input 
to the program is the coordinates of the two points. The output is the distance 
between the two points. Test your program on points for which you know the 
answer. Use your program to fill out the table below. 
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FERMAT’S PRINCIPLE 


Now we are ready to start exploring the idea in geometrical optics. It can be 
stated as follows: 





If a ray of light passes between two points, out of all the 
possible paths that could be taken, the light ray will follow 
the one for which the transit time is a minimum compared 
to nearby paths. 


This is a remarkable statement! How does nature explore all the possible paths 
and direct the light ray along that very particular one such that it will pass be- 
tween the two points in the least possible time? This becomes even more remark- 
able when we see that the path selection must be done in advance. Richard 
Feynman, the world renowned physicist, stated that nature seems to be able to 
“smell” the correct path in advance. We cannot explain the how or why of this, 
but must remain content with the fact that it does happen. 


To proceed, we must find how to compute the transit time for a ray of light 
following a particular straight path. If v is the speed of light and d is the distance, 
then t, the transit time is given by 


t=d/v 


The value of v depends upon the medium through which the light is passing. If 
we divide the velocity of light in a vacuum (c = 3 X 10!° cm/sec) by the velocity 
in a medium, the result is 


n = c/v = index of refraction for the medium 


For every medium, we can find the value of n. As examples, water has n = 1.33 
while most types of glass have an index of refraction of about 1.50. 


Now we can return to the transit time, and express it in a slightly different form. 
Since 1/v = n/c we have t = nd/c. If along the actual light path t is a minimum, 
then nd/c must also be a minimum. But since c is a constant, nd must be a mini- 
mum when compared to any other path. The quantity nd is called the optical 
path. If a ray passes through a distance d, in a medium with index of refraction 





Fermat’s Principle 


(1) 


(2) 
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ny, and a distance dy in a medium with index of refraction nj, and so on, the total 
optical path is the sum of the parts, or 


total optical path =n,d, +d,d, +n3;d3; +.. 


The point that we have been driving at is that having the least transit time for a 
ray passing between two points is equivalent to saying that the total optical path 
has the least value. Now we will rephrase the idea of geometrical optics in the form 
most suitable for our discovery exploration. 


A ray of light will follow that path between two points for 
which the total optical path is the smallest as compared to 
nearby paths. 


This is known as Fermat’s principle and was discovered in 1650. We have not 
stated the principle completely above. However, to attempt to be complete at 
this point might hide what we are trying to accomplish. Therefore, we will return 
later and complete the principle. 


REFLECTION 


Consider the diagram in Figure 4. Suppose a ray of light is to pass from point A 
to point B. Quite obviously the least transit time path is on a straight line joining 
the two points. However, suppose we require that the ray pass from A to B by 
reflection from the x axis. If we consider only integer values of x falling between 
the two points, all the possible reflection paths are shown. (We have not considered 
reflection points outside the interval x = 1 to x = 11. However, this is certainly a 
reasonable assumption.) 


14 


(3) 





B(11,6) 





Figure 4 


Exercise 3 — Rays Between Fixed Points 


Write a program which will check each of the possible paths in Figure 4, and find 
the one which would be followed according to Fermat’s principle. Have the pro- 
gram output the value of x where the reflection will take place. Assume that n = 1. 


GEOMETRICAL OPTICS 
Reflection 
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Exercise 4 — Rays Between Changing Points 


Modify the program in Exercise 3 to permit the input of the coordinates of the 
two points. We still want as output the value of x where the reflection will take 
place. Check only integer values of x between the two points. Use only positive 
values of y. Assume n = I. 


In Exercises 3 and 4 you found a way to locate the integer value of x closest to 
the reflection point. Suppose that for a given problem the output is that the 
reflection takes place at x = 3. We have checked only integer values of x and so 
know that the actual value probably isn’t 3. 


However, we also can see that the actual value must fall somewhere between 2 
and 4. If we looked at the values of x of 2.0, 2.1,... up to 3.9, 4.0 we could 
narrow the search down to a zone 0.2 wide in which the true value lies. Then we 
could repeat the process and find the reflection point within 0.02 and so on. 


16 
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Exercise 5 — Precision Computation 


Rewrite your program of Exercise 4 to find the reflection point to M significant 
digits. As input, the program should require the coordinates of the two points, 
and M. Try out the program on a set of points with M = 5. Assume that n = 1. 


Now that we can locate the reflection point as accurately as we please, we can 
proceed to learn more about the reflection process. If at the point of reflection 
we construct a line perpendicular (normal) to the reflecting surface, we can com- 
pute the angle from the normal to the ray, both before and after reflection. If 
X, is the x coordinate of the reflection point, convince yourself that the tangent 
of the angle the ray makes with the normal before reflection is 


tan (0,) = (x; — x1)/y}. (4) 


By the same line of reasoning, we can compute the tangent of the angle the ray 
makes with the normal after reflection. 


tan (02) = (x2 — x,)/y2. (5) 
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Exercise 6 — Angles of Incidence and Reflection 


Write a program to accept as input the coordinates of two points. The program 
should compute and output the angle of incidence (the angle between the ray and 
normal before reflection), and the angle of reflection (the angle between the ray 
and the normal after reflection). Arrange your program to compute to six 
significant figures. Restrict values of y to positive values and as in the previous 
examples, assume n = 1. 


Exercise 7 — Discovery 


Use the program from Exercise 6 to investigate a large number of different prob- 
lems. Print the results in tabular form. 
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Now we are ready to “‘discover’’ the first law which follows from the idea as 
stated in Fermat’s principle. You are probably far ahead by now and know what 
we are driving for. 


Exercise 8 — Generalization 


Generalize the results of Exercise 7 into the form of a law. This law should be 
titled the Law of Specular Reflection. Look up the law in your physics textbook 
and compare to your version. 


Law of Specular Reflection (your generalization). 


Law of Specular Reflection (from textbook). 
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Exercise 9 — Corner Reflection 


Consider the diagram in Figure 5. We want to investigate the path taken if a ray 
leaves point I, is reflected from the line OB, is reflected again from the line OA, 
and finally passes through point 2. Write a computer program which uses Fermat's 
principle to find the exact path which will be followed. Input several sets of 
points and examine the results. See if you can draw a conclusion from your 


results. 


i 
@ (x1, ¥1) 


Z 
@ (x2, y2) 





Figure 5 
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REFRACTION 


In all the exercises to this point we have assumed that n (the index of refraction) 
was a constant along all possible paths of the light ray. Clearly, if the ray passes 
from one medium to another while traveling between our two points, this is no 
longer true. Now, when we use Fermat’s Principle, the optical path will be deter- 
mined both by distances and the values of the index of refraction in the various 
path segments. 


Study the drawing in Figure 6. Here we have shown possible paths from point A 
in air (n = 1) to B in water (n = |.33). The dividing line between the two sub- 
stances is the x axis. We can use one of the programs already developed, with 
minor modifications, to investigate how rays pass from one medium to another. 
Before proceeding, we should note that in general, the index of refraction for a 
substance depends upon the wave length of the light ray. To simplify our dis- 
cussion we will assume that the light rays have a constant wave length, and all 
values of n which are used correspond to that wavelength. 





Figure 6 
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Exercise 10 — Refraction With Precision Computation 


Modify the program in Exercise 5 to investigate the problem in Figure 6. Let M 
(the number of significant digits) = 5. Use as input the coordinates of A and B. 
(Note that the y coordinate of B will be negative.) The program should output 
the x coordinate where the ray passes from air into water as determined by 
Fermat’s Principle. Run the program for several sets of points A and B, and 
record the results. Try to generalize these results. 


Exercise 11 — Refraction Between Different Media 


Run the program in Exercise 10 for each of the cases in the table below. In each 
case, assume that A is at (-10, 10) and B is at (10, -10). Examine carefully and 
try to generalize your results. | 


Medium 2 


Oil 
Water 
Air 
Glass 
Air 
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Exercise 12 — Angles of Incidence and Refraction 


Modify the program in Exercise 11 to print out additionally the angles the ray 
makes with the normal to the boundary between the two media both before and 
after passing from one medium to the next. Run your program for the same input 
data as in Exercise 11. What can you see in the results? 


We are looking for some general relationship which describes how a ray will pass 
from one medium into another. Clearly, the results of Exercise 12 do not show 
any such relationship. However, the index of refraction must have something to 
do with the desired result, as well as the angles. At this point, left to your own 
devices, and with lots of time to engage in trial and error investigations, you could 
probably discover the correct relationship. To shorten the time required, a certain 
amount of direction will be provided. 


Exercise 13 — Discovery 


Write a program which uses the results of Exercise 12 to generate the following 
information for each case. 


n,)@,; n,tané, nN, Ccosé ; n,sin@ ; 


nod nytané> N2 cosé > nj siné > 


- 


In each case, 6, refers to the angle made with the normal prior to the passage 


from medium I to 2, and 0) is the angle after. Examine the results for the relation. 
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Exercise 14 — Generalization 


Put the results of Exercise 13 into the form of a law. The law should be titled 
Snell’s Law. Look up this law in your physics text and compare to your version. 


Snell’s Law (your generalization). 


Snell’s Law (from textbook). 


Exercise 15 — Double Refraction 


Write a program to investigate how a ray of light will pass from point 1 to point 2 
in the diagram below. 


@ (x), Yq) 





@ (x2, ¥2) 
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Exercise 16 — A Complex Problem 


Write a program to investigate how a ray of light will pass from point 1 to point 2 
in the diagram below. Assume that the ray must be reflected from the line y = 0 
before passing from the medium to the left of the y axis to the medium on the 
right. 


(x1, ¥1) (x2, y2) 





Now that the laws of reflection and refraction have been discovered, we should 
return to a complete discussion of Fermat’s Principle as was promised. We will 
now restate the principle in its complete form. 


A ray of light will follow that path between two points for 
which the total optical path is either a maximum, a minimum, 
ora stationary value when compared to nearby paths. 
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We have discussed only cases for which the optical path was a minimum. To 
complete the discussion we should illustrate cases where the optical path is either 
4 maximum or has a stationary value. Suppose we have an ellipse (your instructor 
can show you how to draw an ellipse on the chalkboard and locate its focus 
points) and require that a ray of light leave one of the focus points and then pass 
through the other focus after one reflection from the ellipse. The nature of the 
ellipse is such that the total distance from a focus to any point on the ellipse, 
then to the other focus is a constant. Moreover, at the point at which the reflec- 
tion takes place, the angle of incidence equals the angle of reflection in agreement 
with the law we have just discovered. This is an example of the optical path having 
a stationary value when compared to nearby paths. 


Consider the diagram in Figure 7. 


he 


Figure 7 





Using calculus we can show that if the ray passes from A to C then to B, that the 


optical path will be a maximum if 
+ 
a >R(E 2 =) (6) 
n2 — Ny 


————————— 
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RAY TRACING 


We mentioned earlier that most of geometrical optics is founded upon the twin 
ideas of reflection and refraction. Having discovered the basis for these ideas in 
Fermat’s principle we can proceed to some applications which are particularly 
interesting using the computer. Up to this point we have considered in detail 
only the case where light is either reflected or refracted at plane surfaces. Now 
we will examine what takes place if the surfaces are curved. 





Figure 8 


In Figure 8, we have portrayed a light ray moving from left to right, initially 
making an angle 6, with the horizontal. It then encounters the first curved glass 
surface. If we construct the normal N,; we can relate the angles of incidence and 
refraction and the two indices of refraction by Snell’s Law. 


n,Sin@, = n.Sin@; 


Ray Tracing 


(7) 
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Then the ray passes to the second curved surface where similarly we can write: 


ny Sin@ = n, Sind 


The net effect of the ray passing through the two surfaces is that it has been devi- 
ated from its original angle by ai angle 64. Of course, what we are seeing is the 
action of a lens. We will not pursue this particular question further here as it is 
covered completely in all introductory physics texts. 


We are interested, however, in a special case called the thin lens approximation. 
Here, the curvature of the glass surfaces and the dimensions of the lens are such 
that the lens can be thought to exist in a plane. This leads to the situation por- 
trayed in Figure 9. 





Figure 9 


As before, the ray approuches from the left, passes through the lens, and is deviated 
by an amount 6g. In this case we are interested only in the net effect of the lens, 
i.e., that the incident ray has been deviated. We are not concerned with how this 
happens, what the angles inside the lens are etc. If the ray is deviated towards 
the horizontal (which we will henceforth refer to as the optical axis) it is a con- 
verging lens. If the ray is deviated away from the optical axis, the lens is diverging. 


Suppose now that we have a ray initially parallel to the optical axis which then 
passes through either a converging or diverging lens as shown in Figure 1 0. 


28 


(8) 





Converging Lens Diverging Lens 





(a) (b) 


Figure 10 


In both cases, the distance from the lens to where the ray either crosses or appears 
to have crossed the optical axis is known as the focal length f. The focal length 
for the converging lens is positive, and is negative for the diverging lens. You can 
easily see that as the focal length decreases, the deviation angle increases. 


Now we can write some simple equations which when used with the computer 
become very powerful. First, we will restrict all light rays to those that make small 
angles with the optical axis. These are known as paraxial rays. Next, we will con- 
sider only optical systems made up out of thin lenses and spaces between them. 
There are no restrictions on the number of lenses or spaces. 


The key to ray tracing is the fact that at any given point in an optical system (a 
combination of lenses and spaces) a light ray can be described by two numbers. 
The first is the distance away from the optical axis. The second is the angle the 
ray makes with the optical axis. We will use y to describe the distance, and 6 for 
the angle (assumed to be in radians). 


We must be able to describe the effect of either a lens or a space upon these two 
numbers. The space is easier to describe so we will dispose of it first. Figure 11 
portrays a light ray passing through a distance s along the optical axis. Clearly the 
angle has not changed. 


Figure 11 


GEOMETRICAL OPTICS 
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Thus 6) = 6,. As long as @ is quite small we can say that y2 — y, = s@,. Fora 
space s then, we have the two equations 


62 = 6, and y2 = y; + 841. (9) 


Before we can write the corresponding equations for a lens we must establish 
the relation between the focal length of a lens and the deviation angle 0g. Figure 
12 illustrates the relation we are looking for. 


Lens 





Figure 12 


We know that if the ray was initially parallel to the optical axis, it will cross the 
axis a distance f after the lens (in the case of a converging lens) where f is the 
focal length. If y was the distance of the ray from the optical axis before it en- 
countered the lens, and provided that @g is small we can write 


04 = y/f (10) 


Figure 13 shows the effect of a lens upon an arbitrary ray. 





Figure 13 
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Ray Tracing 
In this case the distance from the optical axis is unchanged, or 
Yo: Ya: 
From the diagram it is clear that 
OQ, =0, - 64. 
But 
O35 =yi/f. 
Therefore 
GO, Sy it. 
Now we can summarize the equations we have been developing. For a space s 
0, =0,, 
(11) 
Yo =i +80. 
For a lens with focal length f 
O5=0, ~ yi/f, 
(12) 
Yo ~ 3¥i1- 
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Exercise 17 — A Single Lens 
Use the equations just developed to follow a light ray which initially has y = I, 


6 =—-0.1 through a distance s = 20 and a thin lens with f = 10. Do the calculation 
by hand and draw a diagram of the results. 


Exercise 18 — Two Lenses 


Repeat Exercise 17, except add a second space, s = 10, and a second lens with 
f = -20. 


32 





GEOMETRICAL OPTICS 


—_—_—_ Eh Sse shader 


Exericse 19 — A Program For N Lenses 


Write a computer program which will carry out the ray tracing done in Exercises 
17 and 18. Arrange your program so that an arbitrary number of spaces and lenses 
can be accommodated. The program should accept the initial values of y and 6, 
and N pairs of numbers in the form (n, A). N is the number of optical elements 
(either spaces or lenses). For a space, haven = 1 and A =s. Fora lense have n = 2 
and A = f. In your program, you can test the value of n at each step to determine 
which set of equations should be utilized. The program should output the values 
of y and 6 at the end of the optical system. Try out your program on the ray: 
y =0,6 = 0.1, N = 4; (1,10), (2,20), (1,20), 2,10). Plot a diagram of the results. 


Exercise 20 — Intersection of Rays 


Trace the rays y = 1, 6 = 0, and y = 1, 6 =—-0.1 through the optical system in 
Exercise 19. Where do the two output rays either intersect, or appear to intersect? 


We can now discuss the question of optical objects and images. An object is a 
source of light rays which are given off in all directions. If two different rays 
which start from the same point on the object either intersect or appear to have 
intersected after passing through an optical system an image of the point on the 
object is formed. If the rays actually intersect at the image point, the image is 
real. If the rays only appear to have intersected, the image is virtual. If the image 
point is on the same side of the optical axis as the corresponding object point, 
the image is erect. Otherwise, the image is inverted. The image can be magnified 
or decreased in size. We will not follow this topic any further as the subject is 
covered thoroughly in introductory physics texts. 


Ray Tracing 
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Exercise 21 — Description of Image 


Let us use an arrow for an optical object. The object is located along a line per- 
pendicular to the optical axis. The point of the arrow is one unit away from the 
optical axis, and the tail is on the axis. An optical system is described by N= 2, 
(1, 20), (2, 10). Use the program in Exercise 19 to trace two different rays origi- 
nating on the point of the object through the system. Describe the resultant 
image. Is the image real or virtual, erect or inverted, magnified or decreased in 
size? Where is the image? 


It will now be wise to recall Fermat’s principle as applied to the optical systems 
we have been studying. No matter how complicated or involved the optical system, 
the path actually taken by a light ray through the system is such that the optical 
path is a minimum (generally). 


Now that we have developed computer tools to handle many interesting problems 
in geometrical optics, your instructor may assign additional problems which you 
either can work with programs you have already written, or with programs you 
are now capable of writing. 
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The final material in this unit should be considered optional in nature. It should 
prove very interesting to those students who have been introduced to the matrix 
notation and fundamental matrix arithmetic. If you have not studied this part of 
mathematics your instructor can suggest references to assist you. 


In matrix notation, the equations for a space take the form 


05 1 O] 6, 
y2 s lljy 


The equations for the lens put into the same form are 


05 1 -1/f}} 6, 
y2 O Illy 


Thus, the space matrix is 


10 
S = 
s ] 
and the lens matrix is 
] -1/f 
ie = 
0 ] 


Now, consider the optical system in Figure 14. Given a ray at position 1, we want 
to compute the description of the same ray when it emerges at position 5. Using 
the space and lens matrices we have the following. 
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Lens 1 Lens 2 


Object 





Figure 14 
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However, the product of four 2 X 2 matrices in the last expression is itself a 2 X 2 
matrix which we can call the system matrix M. Now the whole problem has the 
simple form 


Mout Gin 


=M 
Yout Yin 


where M is the product of all the space and lens matrices in the system, the pro- 
duct taken in reverse order. 


Exercise 22 — A Matrix Program 


Write a computer program using MAT commands to compute the system matrix 
for a series of N components, each component being a space ora lens. Test your 
program on the system in Exercise 19. You should, of course, get identical results. 


It is interesting to note that a complete description of the system is contained in 
the four numbers in the system matrix M. All the information concerning the 
image corresponding to a given object can be determined from these four numbers. 
This concept of a system matrix turns out to be an extremely powerful idea in 
several branches of physics. 


Exercise 23 — Description of Image 


Develop a set of rules using the four numbers in a system matrix to determine if 
the image is real or virtual, erect or inverted, and the magnification. 
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NOTES 
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APPLICATION PROGRAM 


We have included this completely documented Application Program to illustrate 
one manner in which you might use the computer to solve a problem. 





The program is listed with a definition of the problem, 2 flow chart, a line by line 
description of the program, a listing and a run. 


As an application program to illustrate the use of flow charting we will examine a 
harmless looking problem which nevertheless requires a computer for solution. 
The problem is to find the root of the equation 


x = Cos(x) 


that lies between x = 0 and x = 1, correct to N significant digits. 


If we plot y = x, and y = cos(x), then find the x value where they intersect, we 
will have the desired solution. We know that at x = O the curve y = Cos(x) is 
above the curve y = x. Our strategy will be to step across towards x = 1 until the 
two curves cross, then back up one step and decrease the step size by a factor of 
10, and start the process again. We will go through the process N times at which 
point we will know the root to the desired degree of accuracy. 
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INPUT 
VARIABLES 
D,|, X, AND N 





YES 


40 


PRINT 
ANSWER 
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Application Program 





Line Description 


uve | States object of problem and establishes input feature. 


110 

120 Sets variables D, I, X and N. (D is an increment value on X, I generates 
130 + the program N times working toward the intersection of the equations, 
140 consequently decreasing D and producing the desired significance. ) 

150 

151 , 

160 Examines IF statement and computes new values of X 

170 

i Alters D and I values. 

200 

501 Examines IF statement with respect to I and N and goes back through 
10 loop or prints answer. 


999 End. 


41 








PHYSICS 


Hewlett-Packard Computer Curriculum 





100 REM ROOT OF X = COSCX) 
116 INPUT N 

120 LET D=e1 

130 LET I=@ 

146 LET xX=@ 

15@ IF X<COSCX) THEN 168 
151 GOTO 178 

162 LET X=X+D 

161 GOTO 15¢ 

17@ LET X=X-D 

186 LET D=D/16 

196 LET [=I+1 

26@ IF I >= N THEN 214 
201 GOTO 158 

219 PRINT “KOOT AT "“INTCA* 1OtNIZTOTN 
999 END 


RUN 


2 
ROOT AT 073 


KUN 


2 5 
RIOT AT © 139088 
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